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Stabilization of helical macromolecular phases by confined bending 
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By means of extensive replica-exchange simulations of generic coarse-grained models for helical 
polymers, we systematically investigate the structural transitions into all possible helical phases for 
flexible and semiflexible elastic polymers with self-interaction under the influence of torsion barriers. 
The competing interactions lead to a variety of conformational phases including disordered helical 
arrangements, single helices, and ordered, tertiary helix bundles. Most remarkably, we find that 
a bending restraint entails a clear separation and stabilization of the helical phases. This aids in 
understanding why semiflexible polymers such as double-stranded DNA tend to form pronounced 
helical structures and proteins often exhibit an abundance of helical structures, such as helix bundles, 
within their tertiary structure. 

PACS numbers: 64.70.-p, 82.35.Lr, 83.10.Tv, 87.15.Cc 


Helical segments are ubiquitous secondary structures 
occurring in most macromolecular systems. The forma¬ 
tion of helical structures is typically attributed to the 
formation of hydrogen bonds along the backbone of lin¬ 
ear polymers, but it is also known that helices are among 
the few generic geometries that a linelike topology can 
form if an orderingprinciple (such as a many-body con¬ 
straint) is present [l|-l3|. 

In seminal works, Zimm and Bragg (ZB) ii showed 
that the crossover between disordered random coil struc¬ 
tures and ordered helical conformations can be described 
by a one-dimensional Ising-like model. Therefore, while 
short-range cooperativity can lead to structural ordering, 
in the ZB model this process is not a phase transition in 
the strict thermodynamic sense [6|, [7|. However, since 
biologically relevant macromolecules are finite systems 
(on an effectively mesoscopic scale), the thermodynamic 
interpretation of structural transitions in such systems 
must address finiteness effects accordingly [8[. 


Primary effects of cooperativity can be addressed by 
generic effective-potential models that allow for the qual¬ 
itative description of helix-coil transitions [^12|. Dom¬ 
inant nonbonded interactions support the formation of 
tertiary structures including single helice^helix bundles. 


collapsed globules, or random coils |13l42Qj. It has been 
shown recently that the alignment of secondary struc¬ 
tures in a tertiary protein fold can be understood as a 
simple two-state process [21 1. 


For a generic flexible polymer chain, in a crystalliza¬ 
tion process succeeding the chain collapse, ordered struc¬ 
tures emerge that are substantially different from ter¬ 
tiary structures known from realistic biomolecules and 
typically do not possess secondary structures 22|, [HI . If 


bending restraints and nonbonded interaction compete 
with each other, as it is the case in self-interacting semi¬ 
flexible polymers, the ordered structures are known to be 


rodlike bundles or toroids 0 . 

However, less is known about the influence of effec¬ 
tive bending restraints upon transition pathways toward 
helical structures. A systematic analysis of the forma¬ 
tion and separation of helical phases in phase space in 
the presence or absence of bending restraints has not yet 
been performed. In this study, we investigate the rel¬ 
evance of this restraint for the separation of structural 
phases by means of replica-exchange Monte Carlo com¬ 
puter simulations for coarse-grained flexible and semiflex¬ 
ible polymer models. By scanning the spaces of torsion 
parameter strength and temperature, we construct the 
hyperphase diagrams for entire classes of helical macro¬ 
molecules, which allows us to distinguish the different 
pathways to the helical folds and enables us to judge the 
significance of bending restraints in biomacromolecules. 

We employ a generic coarse-grained bead-spring model 
for elastic, self-interacting polymers with torsional inter¬ 
action. The polymer is represented by a linear chain 
of N monomers. The bending energy of flexible poly¬ 
mers is zero. For semiflexible polymers, excitations of 
the bond angle formed by successive bonds are subject 
to an energetic penalty. Torsion is induced by an out- 
of-plane torsion angle between three successive bonds. 
Nonbonded monomers interact via long-range attractive 
and short-range repulsive van der Waals forces, mod¬ 
eled by the Lennard-Jones (LJ) potential. The energy 
of a conformation X = {xi, X 2 ,..., xn}, where x^ is 
the location of the ith monomer, is given in units of 
the LJ energy scale e by E{X)/e = '^Lj(^ 2 j) + 

Hi ^bond('^i Sq H/- '^bend(^/c) The 

dimensionless Lennard-Jones potential with cut-off is 
given by ULj(r) = 4[(cr/r)^^ - (cr/r)^] - Uc if r < Tc = 
2.5cr, where r is the distance between two nonbonded 
monomers and Vc ~ —0.0163 is a constant shift to en¬ 
sure ULj(rc) = 0. For r > Tc, we set ULj(r) = 0. 
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FIG. 1. Scatter plots of conformations of the (a)-(d) 
semiflexible (bending restrained) and (e)-(h) flexible (bend¬ 
ing unrestrained) polymers with 40 monomers in (^ 1 ,^ 2 ) 
space. Light gray regions represent the generalized ensemble 
of all conformations found at all temperatures T and torsion 
strengths St simulated. Black regions correspond to the most 
populated states at given St values. Representative confor¬ 
mations are shown. 

Distances are measured in units of the length scale ro, 
given by the location of the LJ potential minimum. The 
van der Waals radius of a monomer is chosen to be 
(j = 2“^/^ro. The FENE (finitely extensible nonlin¬ 
ear elastic) bond potential is employed in the form [25[ 
'^bond(^) = — [(^ ~ The bond strength 

is Sr = —KBj^ 12 e] the parameters were set to standard 


values K = (98/5)er^ and R = (3/7)ro. The pay-off for 
bending the chain is 'i;bend(^) = 1 —cos(6> —6>o), where Oq is 
the bond angle in the ground state. The bending energy 
scales with sq = 80 / 6 . Eor the simulations of the flexible 
polymer Se = 0 (no bending restraint), whereas for the 
semiflexible polymer Se = 200e was chosen. Eventually, 
the torsion potential is 'r’tor(^) = 1 — cos(r — tq), with 
the dihedral torsion angle r and its equilibrium value tq . 
The relative energy scale is St = 5'r/e. The choice of 
reference angles tq = 0.873 and Oq = 1.742 allows for 
helical segments in the ground-state structures that re¬ 
semble right-handed a helices with about 4 monomers 
per turn. In the following, e, ro, and /cp are set to unity. 
Eor the simulation of polymers with up to 60 monomers, 
replica-exchange Monte Carlo parallel tempering simula¬ 
tions were performed [2614^. 

The propensities of polymers with {Se > 0) and with¬ 
out bending restraint {Se = 0) to form stable helical 
structures are investigated under thermal conditions con¬ 
trolled by the canonical heat-bath temperature T. The 
variation of the torsion strength St enables the study of 
an entire class of helical polymers. Representations of 
transition channels in generalized ensembles have turned 
out to be beneficial [1, Therefore, we discuss 

the folding channels of the helical polymers in the mul¬ 
tiplicative canonical ensemble provided by the parallel 
tempering method. We introduce a pair of order param¬ 
eters which are effectively defined by the average total 
energies per monomer of the nonbonded LJ interactions 
between all monomers and their neighbors up to 6 bonds 
away, 

^ N-2 N 

9i(X)=e—^ E 06,i-iVLj(ry), (1) 

i=l j=i+2 

and all others, 

N-2 N 

92 (X) = e—^ ^ 0i-i.7VLj(rii), (2) 

i=l j=i-\-2 

where 0/^/ = 1 if /c > / and zero otherwise. In a single 
long helix, qi is minimal and q 2 maximal, whereas for 
helix bundles with increasing number of segments q 2 gets 
smaller and qi larger. 

Eor a 40mer, Eig. [1] depicts for a selection of tor¬ 
sion strength values St the distributions of conformations 
found in the generalized ensemble that covers the tem¬ 
perature interval T G [0.1, 2.0]e//cB in (^ 1 ,^ 2 ) space. The 
left column of figures [(a)-(d)] shows the helical tran¬ 
sition pathways for the bending-restrained semiflexible 
polymer {Se = 200e) and the right column [(e)-(h)] for 
the unrestrained flexible polymer {Se = 0). The lightgray 
region represents the area in (gi, ^ 2 ) space, in which con¬ 
formations were found at all temperatures and torsion 
strengths in the simulations. This distribution, which is 
independent of T and 5'^-, gives a first impression of the 
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FIG. 2. Ratio of order parameters {q 2 )/{qi) for lowest-energy 
structures at various values of Sr in the cases of restrained 
and unrestrained bending. 

differences of the conformational phases of entire classes 
of semiflexible and flexible polymers. It depicts the pos¬ 
sible folding channels for the polymers. From the fig¬ 
ure, it is obvious that the distributions spread out much 
more for the bending-restrained polymer. Individual sec¬ 
tions (phases) are clearly separated, with less-populated 
regions in-between. This is different in the case of flex¬ 
ible polymers. Although conformational phases can be 
identified as well, their separation is much less promi¬ 
nent. These differences can be interpreted in the way that 
in the case of semiflexible helical polymers, structural 
phases are more stable, because these are surrounded by 
entropically suppressed regions, which cause free-energy 
barriers and phase separation between the helical phases. 
The black regions in the figures represent the populations 
in (^ 1 ,^ 2 ) for fixed Sr values, i.e., for individual polymer 
systems and confirm that semiflexible polymers with a 
certain torsion strength prefer to form tertiary structures 
inside a distinct helical phase only, which is not the case 
for flexible polymers. 

The differences in their structural behavior are also 
clearly visible when plotting the order parameter ratio 
{Q 2 ) / (^1) for the lowest-energy structures found in the 
simulations, as shown in Fig. [2l The step-like decrease 
of this ratio for the bending-restrained polymers enables 
the location of the threshold values of Sr where structural 
phases are separated. 

For a more systematic analysis of the folding behav¬ 
ior and its dependence on the torsion strength Sr and 
temperature T, we define the free-energy in order pa- 
rameter space by ^ 2 ) = -ksTlog 

where ^ 2 ) = - 9i(X))(5(g^ - 

g' 2 (X)) exp[—F^(X)// cbT] is the restricted partition func¬ 
tion in the space of all polymer conformations {X}. Fix¬ 
ing Sr and T, the free energy Fs^^t{qi^Q 2 ) possesses a 
global minimum at order parameter values ^ q^^^)• 
The ensemble of all conformations with these order pa¬ 
rameter values represents a dominant macrostate of the 


system. The space of macrostates that share a charac¬ 
teristic structural feature such as the number of helical 
segments in a helix bundle forms a structural phase. 

Transitions temperatures for the various model param¬ 
eter settings were identified by standard canonical anal¬ 
yses of extremal fluctuations of energy (specific heat), 
structural quantities (e.g., radius of gyration), and order 
parameters ((^i), (^ 2 ), number of helices). This will be 
discussed in more detail elsewhere [IH . Based on this in¬ 
formation, the (^ 1 ,^ 2 ) space can be separated into regions 
(“structural phases”) as shown in Fig. [3] for bending- 
restrained semiflexible (left figure) and unrestrained, flex¬ 
ible polymers (right figure). Black lines represent folding 
trajectories for several single polymers with fixed torsion 
strengths in the interval Sr G [0, 30] in (gi, ^^ 2 ) space upon 
cooling. All trajectories begin in the high-temperature, 
random-coil phase (upper right corner in Fig.[3l i.e., large 
qi , q 2 values) and propagate toward a helical state by de¬ 
creasing the temperature. The folding channels at given 
Sr values effectively connect free-energy minima (dots) 
at various temperatures. 

The structural phases of flexible polymers are less well 
separated, in which case folding channels, after passing 
a liquid phase, end in the solid amorphous phase. This 
general transition behavior is virtually independent of 
the torsion strength Sr- However, the influence of the 
value of Sr upon helix and helix-bundle formation is sig¬ 
nificant for the bending-restrained, semiflexible polymer. 
For = 0, the behavior is similar to that of the flexible 
polymer, but for torsion strength = 1, it crystallizes 
initially into a three-helix bundle and then undergoes a 
solid-solid transition. It emerges from it as a four-helix 
bundle, which is energetically slightly more favorable at 
very low temperatures. Three-helix bundles clearly form 
at sufficiently large torsion strength (e.g., for Sr = 5). 

Increasing the torsion strength further favors the ex¬ 
tension of helical segments, compared to local nonbonded 
contacts. The torsional interaction over compensates 
what had been an energetic gain of nonbonded monomer- 
monomer contacts despite necessary bending penalties. 
The number of turns is reduced to a single one and a 
double-helix forms in the solid phase {Sr = 8). Bending- 
restrained polymers with Sr = 20 coexist in a transition 
state between single- and double helix, i.e., in an inter¬ 
mediate ordered solid phase a single helix is formed first, 
which upon further cooling splits into a double-helix at 
the expense of the formation of a single turn. This torsion 
strength marks the threshold at which distant monomer- 
monomer contacts are still formed. For torsion strengths 
close to Sr = 30 and beyond, only stable single-helix 
phases form. 

The complete structural hyper-phase diagram, 
parametrized by temperature T and torsion strength 
5 't-, is depicted in Fig. |4] for bending-restrained, semi¬ 
flexible polymers (left) and for unrestrained, flexible 
polymers with torsion (right). The phase diagram 
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FIG. 3. Structural phase diagrams for bending-restrained semiflexible (left) and unrestrained flexible polymers (right) in 
(^ 1 ,^ 2 ) order parameter space for the temperature and torsion strength space (T,Sr) covered in our simulations. Colored 
regions represent structural phases. Black dots locate free energy minima at given T and St values. Trajectories show the 
helical folding pathways at fixed torsion strengths St by decreasing the temperature. 


for the semiflexible polymer exhibits apparently more 
structure in the folded regime at temperatures T <0.5 
over the entire interval of torsion strengths. Whereas at 
torsion strengths St < 7 four-helix bundles, three-helix 
bundles, and amorphous conformations compete and the 
phases sensitively depend on the temperature, two-helix 
bundles and single-helix conformations are clearly 
dominant for St > 7. Remarkably, the liquid (globular) 
phase disappears for sufficiently large torsion strengths 
{St > 15), in which case direct coil-helix transitions 
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FIG. 4. Hyper-phase diagrams of bending-restrained semi¬ 
flexible (left) and unrestrained flexible polymers (right) with 
40 monomers, represented in the space of the torsion strength 
St as a material parameter distinguishing classes of polymers 
and the temperature T as an external control parameter for 
the formation of structural phases. The color code is the same 
as in Fig. O 


occur. Within the range 15 < 5'^- < 27, helix-helix 
(solid-solid) transitions are present, where single helices 
collapse into two-helix bundles by forming a turn. 
Once the torsion strength dominates over nonbonded 
monomer-monomer interaction, i.e., for St > 27, the 
well-known direct transition from random coils into 
single helices occurs. 

Contrarily, the folding process of flexible polymers 
is hardly affected qualitatively by torsional constraints 
[Fig. 0] (right)]. The three phases of random coils, glob¬ 
ular, and amorphous structures are well separate, but a 
helical phase is nonexistent. Furthermore, if bending is 
not restrained, the liquid phase does not disappear and 
thus a helix-coil transition does not occur. 

In this Letter, we have systematically investigated the 
influence of bending restraints upon the formation of sta¬ 
ble helical phases. We determined all structural phases 
for entire classes of flexible and semiflexible polymers 
with torsion. These results were summarized in struc¬ 
tural hyperphase diagrams for both polymer classes. 

The primary result our study is that an effective bend¬ 
ing restraint along the polymer chain is necessary to sta¬ 
bilize helical structures and, in particular, helix bundles. 
Different helical structure types that are separated by 
entropic gaps in conformational space can only be identi¬ 
fied clearly for semiflexible polymers, whereas for flexible 
polymers torsional barriers alone are not sufficient to sta¬ 
bilize individual helical phases. 

The outcome of this study provides evidence for the 
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natural preference and significance of locally ordered he¬ 
lical secondary structures for semiflexible biopolymers, 
which effectively include DNA and most proteins. Our 
results support the understanding of the almost strict 
confinement of bond angles in polypeptides (such as bio¬ 
proteins), which reduces the set of degrees of freedom 
that participate in their functional structure formation to 
dihedral angles. For this reason, it is unlikely that flex¬ 
ible polymers, i.e., polymers without bending restraint, 
can be vital and functional in a biological system. 
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